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I.  INTRODUCTION 

One  problem,  of  two  under  consideration,  is  to  assess  the  performance  of  a  radar  bomb 
director  (RED)  by  estimating  its  Circular  Error  Probable  (CEP)  given  in  units  of  length.1 
based  on  measurements  of  miss  distances  from  the  target  at  (0. 0)  in  actual  bomb  drops 
that  fall  under  an  uncorrelated  bivariate  normal  distribution  (UBND)  with  zero  means  and 
unknown  variances.  Confidence  intervals  (CIs)  for  the  true  CEP  (CEPT),  will  be  needed, 
where  the  CEPT,  by  definition,  is  the  radius,  RT,  of  the  circle  centered  at  (0.0),  which 
contains  50%,  PCEP  =  0.50,  of  the  distribution.  The  computer  program,  CONFREG, 
discussed  in  Section  VI,  allows  arbitrary  PCEP  values  in  (0, 1). 

In  the  second  problem,  it  is  assumed  a  radar  contractor  is  to  be  awarded  a  bonus  if  the 
computed  CEP  is  below  an  agreed  nominal  value  (RTn,  for  example  by  10%),  and  is  to  be 
penalized  if  over  the  RTn  by  a  like  percentage.  Thus,  it  is  important  to  know  the  level  of 
confidence,  under  these  circumstances,  that  can  be  placed  on  the  CEP,  R.,  computed  from 
experimental  data  as  an  estimate  of  the  RT. 

In  May  1971,  Dr.  M.  P.  Jarnagin  of  the  Naval  Weapons  Laboratory  at  Dahlgren,  Virginia, 
issued  an  informal  memorandum  [1]  in  which  he  described  statistical  procedures  for  address¬ 
ing  these  problems.  Because  of  limited  computing  facilities,  his  analysis  used  approximations 
throughout  and,  consequently,  the  work  was  limited  in  scope.  Our  objective  is  to  replace 
the  approximations  with  precise  calculations  and  to  generalize  the  results  to  some  extent. 
In  many  parts  of  this  report,  excerpts  are  taken  directly  from  his  memorandum. 

Other  studies  addressing  the  first  problem  have  been  carried  out  in  [2]  and  [3].  Again, 
their  analyses  are  based  on  approximations  without  error  bounds  rather  than  on  precise 
calculations. 

A  set  of  hypothetical  individual  bomb  drops  is  shown  in  Figure  1. 


Figure  1.  Typical  Bombing  Pattern,  n  =  22 

1In  this  report,  the  unit  of  length  will  be  meters. 
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Lot  n  denote  the  total  number  of  bomb  drops.  The  target  is  a  point  target  at  (0.0),  the 
origin  of  coordinates.  The  flight  path  of  the  aircraft  dropping  each  bomb  is  along  the  x-axis. 
The  y-axis  passes  through  the  target,  is  perpendicular  to  the  x-axis,  and  is  positive  above  the 
x-axis.  The  x  and  y  coordinates  of  the  impact  point  of  a  given  bomb  drop  are  the  components 
of  the  miss  distance  in  range  and  deflection,  respectively,  in  the  usual  terminology. 

In  Figure  1.  n  =  22,  Bomb  No.  21  hits  at  point  P21  with  coordinates  (in  meters)  x  = 
8.668.  y  =  23.298.2  The  estimated  CEP,  or  more  briefly  R,  is  the  radius  of  the  circle  that 
contains  50%  of  the  impact  points.  An  indicated  CEP  in  Figure  1  is  R  =  26.5.  For  such 
finite  samples,  R  is  not  uniquely  defined,  but  for  a  theoretical  bivariate  normal  distribution, 
RT  is  a  precisely  defined  unique  number. 

The  assumption  is  made  that,  the  x-coordinates  of  the  bomb  drops,  {x,}=  Xi, x2,  •  •  •  ,  x„, 
constitute  a  random  sample  from  a  normal  distribution  with  mean  zero  and  with  unknown 
standard  deviation  a> .  Similarly,  it  is  assumed  that  the  y-coordinates  of  the  bomb  drops, 
{yi}=  yi,y2,  •  •  •  ,y»,  form  a  random  sample  from  a  normal  distribution,  independent  of  the 
distribution  of  the  Xj’s,  with  mean  zero  and  unknown  standard  deviation  ay.  Thus,  the 
impact  points  (xj.  Vi)  are  assumed  to  be  from  a  UBND.  The  results  in  an  actual  test  may  not 
indicate  zero  means.  That  would  imply  a  bias  in  the  performance  of  the  RBD  that  would 
penalize  the  contractor.  Subsequently,  we  allow  for  nonzero  means,  which  extends  the  model 
in  [1], 

The  CEPT(=  RT)  is  a  function  of  the  standard  deviations  ox  and  oy.  We  compute  an 
estimate  R  for  RT  on  the  basis  of  the  experimentally  determined  estimates  sx  and  sy,  for  rrx 
and  <7y,  where 

sy  ~  ,  f=n,  (Assume  x  =  0,  y  =  0),  (1) 

Sx  =  \/Er=i(xi-x)2/f,  sy  =  y/5Z"=i(yi  -y)2/f  >  f=n-l,  x/0,  y /O,  (2) 

and  where  f  denotes  the  number  of  degrees  of  freedom  in  the  data.  The  means  in  (2)  are 
given  by  n  n 

x  =  £  Xj/n,  y  =  y‘/n-  (3) 

i=l  i=  1 

The  pertinent  function  connecting  the  CEP  and  the  estimated  standard  deviations  taken 
from  (1)  is  the  Generalized  Circular  Error  Function  (GCE),  V(K,  c),  where  K  =  R./sx,  c  = 
Sy / Sx ,  [4]  —  [6].  If  the  estimated  standard  deviations  are  taken  from  (2),  then  the  applicable 
function  is  the  Elliptic  Coverage  Function  (ELP),  P(R,  H.K,  sx,  sy).  Both  functions  will  be 
discussed  in  the  next  section.4 

2 For  clarity,  numerical  results  will  be  included  with  the  general  subject  matter. 

3 Note,  in  general,  sx  ^  sy. 

4The  arguments  K  in  V(K,c)  and  P(R,  H,K,  sx,  sy)  have  different  meanings. 
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II.  V(K,  c)  AND  P(R,  H.  K,  sx,  sy)  PROBABILITY  FUNCTIONS 

The  GCE,  also  known  as  the  elliptical  normal  function,  defines  a  probability  function 
V(R/u,  v/u).  GCE  gives  the  probability  of  a  shot  falling  in  a  circle  of  the  xy-plane,  of  radius 
R  and  centered  at  the  origin,  under  a  UBND  with  mean  zero  and  standard  deviations  u,  v. 
This  probability  is  given  by 


1  r  f 

V<R-v>=2 


R  rV  R2— x2 


1  . .  X . 


exp{-^[(-)2  +  (-)2]}dydx, 
2  u  v 


(4) 


where  u  =  ctx,v  =  ay.  Transforming  to  polar  coordinates,  with  x  =  rcos#,  y  =  rsintf  yields 


,  .  1  fH  ,  l..rcos0.9  .rsinfl.o..,  , 

V(R,  u,  v)  =  - -  /  /  exp{  — -[(  —  )2  +  (  —  )2]}  r  d#dr. 

Z7TUV  Jo  Jo  Z  U  V 


Using  elementary  trigonometric  identities 

Y  /»R/u  rn 


V(R,  u,  c)  =  —  [  I  exp[—  ^r2(B  —  Acos0)]  r  d#dr, 
nc  Jo  Jo  2 


with 


A  = 


1  —  c2 


2c2  ’ 

where,  without  loss  of  generality,  v  <  u, 


B  = 


1  +  c2 

T^2"’ 


(5) 


(6) 


(7) 


c  =  v/u,  0  <  c  <  1. 


(8) 


Using  the  fact  that 

I0(x)  =  —  f  exp(x  cos  6)  d0,  (9) 

n  Jo 

where  I0(x)  denotes  the  modified  Bessel  function  of  the  first  kind  and  zero  order  [7,  p.  375], 
one  obtains  from  (6)  and  (9) 


1  /'R/u  B  r2  A  r2 

V(R/u,  c)  =  -  /  exp(-— -)I0(— )rdr. 

c  Jo  2  2 

(10) 

Now  letting 

* 

III 

33 

C 

3 

II 

to 

to 

(11) 

(10)  reduces  to 

i  rK2/2 

V(K,  c)  =  -  /  exp(— B  w)  I0(A  w)  dw. 

c  Jo 

(12) 

The  significance  of  the  V(K,c)  function  in  the  present  context  is  brought  out  in  Figure  2. 
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y 


Figure  2.  An  Elliptical  Distribution  With  Zero  Means 


Suppose  that  x  and  y,  the  rectangular  coordinates  of  the  impact  points  of  a  series  of 
weapons,  are  distributed,  each  with  zero  mean,  and  with  standard  deviations  ax  and  ay  with 
oy<ax.  The  ellipse  in  Figure 2  has  these  standard  deviations  for  its  semi-axes  lengths  but 
has  no  further  significance.  Now  consider  a  circle  of  radius  R  =  Kdx  meters,  K  being  a 
dimensionless  ratio.  What  fraction  of  the  weapons  of  this  distribution  hits  in  the  interior  of 
this  circle?  Or  equivalently,  as  noted  earlier,  what  is  the  probability  that  a  single  weapon 
from  the  distribution  hits  inside  the  circle?  The  required  probability  is  given  by  the  function 
V(K,c)  of  (10)  or  (12),  where  K  can  take  all  values  from  zero  to  infinity  and  0  <  c  <  1.  If 
c  >  1,  then  <7X  and  oy  are  interchanged  with  no  lose  of  generality. 

The  problem  considered  here  is  of  an  inverse  character.  Estimates  sx  and  sy,  for  rrx  and 
ct y ,  are  obtained  from  bomb  drops  as  noted  in  (1).  Then  K  is  obtained  from  a  subroutine 
INVGCE,  [8],  and  hence  R(=  Ksx),  such  that  V  =  0.50  and  c  =  sy/sx.  The  resulting  value 
of  R.  the  radius  of  the  circle,  which  contains  50%  of  the  distribution,  is  the  computed  CEP. 
Quick  estimates  for  K  can  be  found  from  the  tables  included  in  [4],  [5],  and  [8]. 

The  analysis  in  [1]  is  extended  here  to  include  nonzero  means  (x  /  0.  y  /  0).  For  this 
purpose,  ELP  is  needed.  ELP  specifies  a  probability  function  P(R,  H,  K.  u,  v),  [6]  and  [9], 
where  u  =  (Tx,v  =  ery,  that  gives  the  probability  of  a  shot  falling,  under  a  UBND  with  mean 
zero  and  standard  deviations  u,  v,  in  a  circle,  T,  with  radius  R  and  centered  at  (H.  K)  of  the 
xy-plane.  This  probability  is  given  by 


P(R 


/•n-ffi  /*l 

:,  H.  K,  u,  v)  =  - -  /  / 

2ttu v./h_r  ./ K 


H+R  /•K+v/R2-(x-H)2 


VR'2-(x-H)2 


exp 


<n)2  +  (v)2 


dydx.  (13) 
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Introducing  dimensionless  variables,  let  x  =  (\/2u)s  and  y  =  (\/2v)t,,  then  (13)  becomes 

i  /*h+r  /»k-f  (u/v)-^/  r2  — (s— h)2 

P(r,  h,  k,  u/v)  =  -  /  exp(-s2)  /  _  exp(-t2)dtds  (14) 

7T  ./h— r  «/k— (u/v)y^r2  —  (s— h)2 

with  r  =  R/(\/2u),  h  =  H/(\/2u),  andk  =  K/(\/2v).  Finally,  (14)  can  be  written  as 

P(r,  h,  k,  u/v)  =  ^-j= 

where 

Fl(s)  =  exp(— s2)  aerf  [k,  (u/v)\/r2  —  (s  —  h)2  ],  (16) 

an<*  aerf  (a,  b)  =  -^=  [  exp(— z2)dz.  (17) 

V7r  J  a— b 

Equation  (15)  is  used  for  computational  purposes  as  described  in  [9],  but  using  simple 
linear  transformations  in  (13)  it  is  easy  to  show  that  ELP  can  also  be  interpreted  as  a 

probability  function  P(R,  H,  K,  u,  v),  which  gives  the  probability  of  a  shot  falling  under  a 

UBND  with  mean  (H,  K)  and  standard  deviations  u,  v,  in  a  circle  of  radius  R  centered  at 
the  origin  of  the  xy-plane.  It  is  in  this  context  when  reference  is  made  to  ELP. 

The  problem  considered  here,  as  above  for  the  GCE  case,  is  of  an  inverse  character. 
Estimates  sx  and  sy,  for  ax  and  <ry,  are  obtained  from  bomb  drops  as  noted  in  (2).  Then  R  is 
obtained  from  the  subroutine  INVELP,  [10],  with  P  =  0.50,  u  =  sx,  and  v  =  sy.  The  resulting 
value  of  R,  the  radius  of  the  circle  in  meters,  which  contains  50%  of  the  distribution,  is  the 
computed  CEP  for  the  case  of  nonzero  means.  Quick  estimates  for  R  can  be  found  from  the 
tables  included  in  [9]  and  [10]. 5 

III.  COMPUTATION  OF  THE  ESTIMATED  CEP 

The  data  determined  from  actual  bomb  drops  and  used  as  inputs  for  computing  an  esti¬ 
mate  of  the  CEPT  are  the  quantities  sx  and  sy,  the  estimates  of  the  true,  but  unknown,  stan¬ 
dard  deviations  ox  and  ay.  If  n  is  the  number  of  bombs  dropped,  and  if  (xi,  yQ,  •  •  •  ,  (xn,  yn) 
are  the  coordinates  of  the  n  observed  impact  points  (see  Figure  1),  it  is  assumed  that  the 
x/s  and  the  y/s  constitute  samples  from  independent  normal  distributions  with  means  zero 
or  nonzero  and  standard  deviations  ax  and  ay. 

No  attempt  is  made  to  estimate  any  correlation  coefficient  that  may  be  indicated  in  the 
data,  although  such  an  inclusion  would  present  no  difficulty.  Also,  although  s2  and  s2  give 
unbiased  estimates  for  cr2  and  cr2,  the  estimates  for  <xx  and  ay  given  by  (1)  and  (2)  are  not 
unbiased;  nevertheless,  they  are  invariably  used. 

Using  (1)  and  (2),  we  obtain  from  the  data  given  below  (see  Figure  1)  values  for  the  CEP 
by  using  the  subroutines  INVGCE  and  INVELP. 

5The  roles  of  r  and  R  are  interchanged  in  [10]  from  their  use  herein. 
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The  made-up  data  was  obtained  from  a  subroutine  that  generated  normally  distributed 
random  numbers  with  zero  mean  and  rrx  =  35  and  then  run  again  with  oy  =  15.  The  data 
generated  is  shown  in  Table  1  and  graphically  in  Figure  1. 

Table  1.  Normally  Distributed  Data  oF  Figure  1.  n  =  22 


i 

Xj 

Vi 

i 

Xi 

.Vi 

i 

34.281 

6.5922 

12 

-5.7568 

23.368 

2 

-25.493 

-5.5825 

13 

22.839 

-16.041 

3 

-15.329 

-31.638 

14 

-3.3694 

-5.9377 

4 

14.245 

4.8120 

15 

34.301 

-5.1021 

5 

-45.153 

.9773 

16 

3.1264 

-7.5556 

6 

-9.2964 

7.4277 

17 

28.403 

.8495 

7 

47.826 

4.5984 

18 

-14.379 

3.2486 

8 

22.682 

-2.8901 

19 

-48.324 

.7353 

9 

33.730 

-16.422 

20 

-41.711 

-19.780 

10 

-19.193 

10.269 

21 

8.6679 

23.298 

11 

-30.543 

-11.832 

22 

14.292 

-18.365 

Table  2  contains  the  quantities  required  to  obtain  the  listed  CEPs  using  the  data  from 
Table  1.  The  first,  two  columns  use  (1)  with  subroutine  INVGCE;  the  last  two  columns  use 
(2)  with  subroutine  INVELP,  where 

c  =  sy/sx,  V  =  P  =  0.50,  H  =  x,  K  =  y.  (18) 

Table  2.  CEPs  Based  on  Table  1  Data 


x  =  0 

o 

II 

l>> 

x  =  0.26571 

y  =  -2.4986 

sx  =  27.592 

sy  =  13.322 

sx  =  28.240 

sy  =  13.393 

c  =  .48281 

CEP=23.717 

c  =  0.47427 

CEP=24.269 

V  =  P  =  PCEP  =  0.50  throughout;  but,  as  mentioned  earlier,  the  subroutine  CONFREG 
of  Section  VI  allows  any  value  in  (0, 1)  for  V  or  P. 

The  next  section  contains  an  analysis  to  find  CIs  for  the  CEPT. 

IV.  CONFIDENCE  INTERVALS 

In  this  section,  we  describe  three  types  of  confidence  regions:  conventional  (CCI),  mini¬ 
mum  (MCI),  and  symmetric  (SCI)  for  the  true  parameter  point  (<rx,  ay)  in  half  of  the  quarter 
plane  of  the  erx-<7y  plane  (0  <  <ry/ax  <  l).6  Using  these  results,  the  three  CIs  (CCI,  MCI, 
and  SCI)  are  established  for  the  true  value  of  the  CEP  (denoted  earlier  by  CEPT  or  RT, 
where  an  estimate  is  denoted  by  CEP  or  R).  Hence,  our  objective  is  to  get  meaningful  CIs 

6  It  is  assumed  «7y<  <7X  without  loss  of  generality. 
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for  the  CEPT.  The  approach  to  achieving  this  objective  is  by  analyzing  the  confidence  that 
can  be  placed  in  the  estimates  sx  and  sy,  which  are  quantities  deduced  directly  from  the  raw 
experimental  data.  This  Jarnagin  idea  forms  the  basis  for  all  results  that  follow. 

Also,  the  results  in  actual  tests  may  not  appear  to  be  compatible  with  zero  means,  but  we 
follow  the  discussion  in  [1]  and  proceed  first  assuming  zero  means,  which  makes  the  analysis 
easier  to  follow.  Nevertheless,  most  of  the  arguments  used  also  apply  to  the  case  of  nonzero 
means  except  for  obvious  deviations  such  as  using  (2)  instead  of  (1)  and  using  the  inverse  of 
ELP,  INVELP,  instead  of  the  inverse  of  GCE,  INVGCE. 

Let  a  confidence  level  a  =  1  —  \/PC,  where  confidence  regions  are  determined  with  an 
assigned  probability,  PC;  say  for  illustrative  purposes,  PC  =  (1  —  a)2  (=  0.90).  Then  a  Cl 
for  ax  with  probability  \/PC  =  (1— a)  (=  0.9486833)  is  obtained  by  finding  appropriate  inte¬ 
gration  limits  for  the  integral  of  the  chi-square  density  function,  the  required  distribution  for 
the  variate  (ns2/er2)  [11,  pp.  276-286],  where  the  integral  is  evaluated  at  a/2  (=  0.02565835) 
and  (1  —  a/2 )  (=  0.97434165).  Mathematically,  we  have 

C (xl/2,  xl-a/2),  f)  =  C(0,  Xa-o/2),  0  -  C(0,  xl/2,  f).  (19) 

where  C(7,  6 ,  f)  denotes  the  chi-square  distribution  function  with  f  degrees  of  freedom, 
namely  ,  rS 

c (7,  f)  =  2f/2~r(F/2)  J  e_Z/2  zf/2_1  dz’  0  ^  O'  <  s  <  °°-  (2°) 

Then  \2a/2  and  x21_q/2),  depending  on  both  a  and  f,  are  determined  by  solving  the  inverse 
problems 

C(0,x«/2,  f)  =  a/2  (=  0.02565836),  C(0,  X(i-a/2)»  f)  =  C1  -<*/2)  (=  0.97434165)  (21) 

and  by  using  the  subroutine  DGINV,7  [6].  For  the  numerical  examples,  using  (21),  one 
obtains 

xl/2  =  11.028377,  xfi—a/2)  ~  36.677187,  using  (1)  with  f  =  n  =  22,  (22) 

Xq/2  =  10.327293,  X(i_a/2)  =  35.377041,  using  (2)  with  f  =  n  -  1  =  21.  (23) 

Introducing  the  notation 

CLx  =  Sxv/f/xji-a^p  CHx  =  s*>/f/Xa/2»  CLy  =  sy^f/X(i-a/2).  CHy  =  Syx/f/y2^,  (24) 
with  (22),  (23),  the  CCIs  for  ox  are  established,  [11,  p.  276], 

21.3698  =  CLx  <  crx  <  CHx  =  38.9711,  use  (1)  with  f  =  n  =  22,  (25) 

21.7579  =  CLx  <ax<  CHx  =  40.2703,  use  (2)  with  f  =  n  -  1  =  21.  (26) 

7DGINV  finds  x  of  the  incomplete  gamma  function,  P(a,x),  [7,  p.260],  [12],  [13],  given  P  and  a,  where 

C(0,  x2,  f)  =  P(a,  x)  =  ^  /0X  e-t  ta~ 1  dt,  x  =  X2/2,  a  =  f/2. 
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Similarly,  CIs  are  obtained  for  ay  using  relations  in  (24),  so  that, 

10.3176  =  CLy  <  <Ty  <  CHy  =  18.8157,  use  (1)  with  f  =  n  =  22.  (27) 

10.3190  =  CLy  <  ay  <  CHy  =  19.0988,  use  (2)  with  f=n  —  1  =  21.  (28) 

We  note  that  a  PC  of  0.9  may  be  fairly  high;  a  lower  level  would  result  in  smaller  CIs  for 
<rx  and  ay. 

It  is  assumed  throughout  that  the  ct’s  for  the  ax  and  oy  CIs  are  equal.  For  the  more  general 
situation,  they  would  not  be  required  to  be  the  same.  In  that  case,  subroutine  CONFREG 
and  a  few  supporting  routines  would  require  a  number  of  changes,  but  the  changes  would  be 
straightforward. 

In  the  earlier  discussion,  the  CCI  with  an  assigned  probability  of  \/PC  was  defined  by 
C(X(i-«/»),  00,  f)  =  C(0,  xl,2-  f)  =  (1  -  v/PC)/2,  C(0,  X(,-o/2),  f)  =  (1  +  VPC)/2.  (29) 

The  CIs  in  this  case  for  <rx  and  oy  are  given  by  (25)  and  (27)  or  (26)  and  (28). 

Referring  to  Figure  3,  the  vertical  stripped  region,  where  (25)  holds  and  ay  extends  from 
zero  to  infinity,  indicates  the  (1  —  a)  confidence  region  for  <rx.  The  horizontal  stripped  region, 
where  (27)  holds  and  crx  extends  from  zero  to  infinity,  indicates  the  (1  —  a)  confidence  region 
of  ay.  The  intersection  of  the  two  regions  generates  a  rectangle,  shown  in  the  figure  as  the 
crosshatched  area,  with  the  lower  left-hand  corner  having  the  coordinates  (CLx,  CLy)  and 
the  upper  right-hand  corner  with  coordinates  (CHx,  CHy). 
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1  1  1  1 
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0  1  1 - 1  *■  1-1  l-i  1  l  -  i  I  i - o 

0  5  10  15  20  25  30  35  40  45* 


Figure 3.  90%  Confidence  Region  For  (<rx,  ay) 

This  intersection  specifies  the  PC  (=  90%)  conventional  confidence  region  for  the  point 
(<tx ,  Oy)  based,  in  general,  on  the  data  of  which  Table  1  is  an  example.  This  follows  because 
it  is  assumed  that  the  x  and  y  coordinates  of  the  bomb  drops  each  arise  from  independent 
normal  distributions  with  unknown  standard  deviations,  <rx,  ay,  respectively.  Hence,  the 
product  of  the  specified  probabilities  (1—  c*)*(l—  a)  (=  0.9486833*0.9486833  =  0.90)  =  PC. 
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The  rectangular  reeion  for  the  zero  means  ease  is  shown  in  Fiorure  4. 


Figure 4.  90%  Confidence  Region:  uses  (1),  n  =  22 
Similar  arguments  result  in  a  (90%)  rectangular  confidence  region  for  ax  and  ay,  using 
(26)  and  (28),  the  nonzero  means  case. 

Noting  that  Xa/2and  X(i_q/2)  are  not  unique  for  a  given  a,  they  can  be  varied  to  obtain, 
from  the  subroutine  CONFMS,  two  other  CIs,  namely  MCI  (a  minimized  Cl)  and  SCI  (a 
symmetrized  Cl).  It  is  understood  that  VT/xf/2  an<^  ^f/xp-a/2)  are  the  same  for  the  <xx-CI 
and  the  corresponding  <ry-CI.  For  MCI,  the  quantity 

CHx  -  CLx  (  see  (24)  )  (30) 

is  minimized,  which  means  CHy  —  CLy  is  also  minimized.  In  the  case  of  SCI,8  the  quantity 

|CHx  +  CLx  -  2|  =  Wilxln  +  'Ji/xl-am  -  21,  (3i) 

is  minimized,  which  implies 

|CHy  +  CLy  -  2|  =  |\Zf/xl/2  +  Vi /x?.-„/2)  -  2|  (32) 

Sy 

is  also  minimized.  It  should  be  noted,  however,  that,  for  any  specific  data  set,  symmetrization 
is  not  possible  for  some  a,  if  %f/Xg/2  >  2,  where  C(0,  X5/2’  0  =  <*/2.9 

The  results  for  ax  and  ay  CIs,  using  the  data  from  Tables  1  and  2,  as  well  as  results  from 
the  discussion  above,  are  summarized  below  in  Table  3. 

8 Note  that  with  (25),  the  CCI  is  far  from  symmetrical;  CLx  is  23%  below  sx,  whereas  CHx  is  41%  above  sx. 

9Since  a  quantity  b  =  \/f/Xa  <  v^f/xj/ 2  is  needed  in  the  next  section,  CONFREG  requires  b  <  2  for  symmetrization. 
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Introducing  notation  similar  to  (24),  but  referring  to  MCI  and  SCI  with  corresponding 
appropriate  values  for  \^2  and  _Q/2) : 

MLx,  SLx  =  sxv'f/x?1_Q/2),  MHx,  SHx  =  sxv/f7S, 

(33) 

MLy,  SLy  =  sy^f/x^_Q/2),  MHy,  SHy  =  svv/f/X“/2  - 
Table 3.  Listing  of  CCI.  MCI,  and  SCI  Using  Data  From  Table! 


x  =  0 

x  =  0.2657 

y  =  0 

y  =  -2.4986 

sx=27.5922 

sx=28.2402 

sy=13.3219 

sy=13.3934 

CLx 

21.3698* 

21.7579* 

CLy 

10.3176* 

10.3190* 

CHx 

38.9711* 

40.2703* 

CHy 

18.8157* 

19.0988* 

MLx 

20.5596 

20.8941 

MLy 

9.9264 

9.9094 

MHx 

37.6151 

38.8057 

MHy 

18.1610 

18.4042 

SLx 

18.3814 

18.5228 

SLy 

8.8748 

8.7847 

SHx 

36.8031 

37.9576 

SHy 

17.7689 

18.0020 

Table  3  from  the  3rd  through  the  8th  rows  and  the  corresponding  six  columns  can  be 
looked  at  as  a  6x0  matrix  with  elements  a(i,j),  where,  for  example,  a(l,l)— >  CLx.  Hence, 
a(l,2)  lists  the  left  endpoint  of  CCI  for  the  x  =  (lease;  the  right  endpoint  is  a(2,2), 
thus  CCI  =  (21.3698.  38.9711).  The  MCI  for  the  y  =  -2.4986  case  is,  from  the  ta¬ 
ble,  (a(3, 6),  a(4, 6))  =  (9.9094,  18.4042).  Asterisked  quantities  were  given  previously  in 
(25)  -  (28). 

We  are  now  in  a  position,  following  [1],  to  determine  the  CIs  for  the  CEPT.  Starting 
with  the  zero  means  case  as  specified  by  (25)  and  (27),  consider  the  <rx  —  oy  quarter  plane 
(cTy  >  0,  <7X  >  0)  shown  in  Figure  4,  and  note  that  nonintersecting  curves  of  constant  R, 
which  are  concave  towards  the  origin,  can  be  drawn.10  These  constant  R-curves,  which  are 
not  circular  segments,  are  symmetric  about  the  45  degree  line,  c  =  1 ,  but  we  need  to  consider 
only  half  of  the  quarter-plane  with  oy/ox.  <  1. 

Then  a  straight  line  L  emanating  from  the  origin  satisfies  the  equation  oy  =  M  ax,  M 
a  constant.  For  the  point  ox  =  sx,  ay  =  sy,  we  have  0  <  M  =  c  =  sy/sx  <  1.  Again, 
see  footnote  10  below:  R-curves  are  definitely  not  segments  of  circles.  Typical  curves  are 
shown  in  Figure  4.  They  are  increasing  in  R.  as  they  move  away  from  the  origin.  Consider 
a  line  L,  then  since  V  =  1/2  and  c  is  fixed,  it  follows  K  is  also  constant  on  L.  Therefore,  as 
ax  increases,  moving  along  L,  R  must  also  increase  on  L  to  retain  K  constant  on  L. 

10These  constant  R-curves,  used  only  for  illustration,  can  be  obtained,  referring  to  the  ax  —  cry  plane,  as  follows:  For  a  fixed 
value  of  R,  determine  a  value  of  K  =  Ki(=  R/<rxi)  from  V(Ki,0)  =  1/2  =  erf(Ki/\/2),  [5,  8, and  7  (p.297)].  This  determines 
the  crxi  such  that  for  a  given  smaller  crx  =  ax 2,  which  fixes  a  K  =  K2,  a  value  of  C2  is  found  by  solving  the  inverse  problem 
V(K2,C2)  =  1/2  that,  in  turn,  determines  crV2  thus  this  is  another  (aX2,  cry2)  point  of  the  constant  R- curve  in  the  <rx  —  ay  plane. 
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Note  that  the  point  (sx,  sy)  is  on  the  middle  R-curve  of  Figure 4,  with  R  =  CEP,  inter¬ 
sected  by  the  line  with  slope  c  =  sy/sx  =  0.483.  In  addition,  the  R-curve  containing  the 
point  (CLx,  CLy),  which  is  also  on  L,  (CLy/CLx  =  sy/sx)  has  an  R-value,  CLR,  given  by 
solving  K  =  CEP/sx  =  23.717/27.5922  =  CLR/CLx.  Thus,  CLR=18.3685,  where  CLx  is 
taken  from  Table  3.  Similarly,  considering  the  R-curve  passing  through  (CHx,  CHy)  on  L, 
with  K  =  CEP/sx  =  CHR/CHx,  and  solving  for  CHR,  one  obtains  CHR=33.4997.  Thus, 
we  take  as  the  CCI  for  CEPT,  CCR  =  (CLR,  CHR)  =  (18.3685,  33.4997).  The  MCR 
and  SCR  for  R  are  obtained  in  the  same  way  and  defined  by  MCR  =  (MLR.  MHR)  and 
SCR  =  (SLR,  SHR).  The  results  are  summarized  in  the  first  and  second  columns  of  Table  4. 

In  the  nonzero  means  case,  using  (26)  and  (28),  H  and  K  are  held  fixed  throughout. 
The  distribution  density  function  is  lower  and  flatter  on  each  L  line  as  ax  and  oy  increase 
with  <Ty/crx  constant;  therefore  R  must  increase.  Also,  R/sx  is  not  constant  along  an  L 
line  since  the  normalized  variables  h,  k  vary  along  L,  and  consequently,  R  / sx  also  varies 
since  P  =  0.5;  see  (14).  Therefore  CLR(CHR)  must  be  obtained  using  the  subroutine 
INVELP,  with  arguments  P  =  1/2,  H  =  x,  K  =  y,  ax  =  CLx(CHx),  and  ay  =  CLy(CHv). 
Hence  the  CCR  =  (CLR,  CHR)  for  the  nonzero  means  case  is  (18.7754,34.5007),  where, 
from  Table  3,  CLx=21.7579,  CLy=10.3190,  CHx=40.2703,  and  CHy=19.0988.  Note  the 
MCR  =  (MLR,  MHR)  for  the  CEPT  in  this  nonzero  case  is  not  precisely  minimized,  nor 
is  SCR  =  (SLR,  SHR)  precisely  symmetrized.  The  overall  numerical  results  are  listed  in 
Table  4. 

Table  4.  Listing  of  CCI,  MCI,  and  SCI  for  CEPT 
Using  Data  From  Tables  1  and  3 


Clf  Means  — ► 

x  =  0,  y  =  0 

x/0,  y/0 

CLR 

18.3685 

18.7754 

CHR 

33.4977 

34.5007 

MLR 

17.6720 

18.0454 

MHR 

32.3322 

33.2535 

SLR 

15.7998 

16.0451 

SHR 

31.6342 

32.5316 

Recalling  the  second  objective  (see  page  1,  2nd  paragraph)  of  establishing  whether  the 
design  contractor  of  the  RBD  is  to  be  rewarded  or  penalized,  two  other  confidence  regions 
are  introduced.  They  form  the  discussion  in  the  next  section  with  the  idea  of  sequencing 
bomb  drops.  We  follow  the  reasoning  in  [1], 
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V.  CONFIDENCE  INTERVALS  BY  SEQUENTIAL  TESTING 

Let  us  assume  we  are  dealing  with  the  zero  means  ease:  (25)  and  (27).  Begin  by  assigning 
RTn,  an  agreed  nominal  value  for  the  CEP.  and  a  fixed-decimal  FP.  Then  assume  that 
if  the  CEP  <  (1  —  FP)RTn  -  RN.  within  the  assigned  confidence  level,  the  contractor  is 
rewarded  (COR)  and  if  CEP  >  (1  +FP)RTn  =  RP.  the  contractor  is  penalized  (COP).  With 
a  numerical  example,  suppose  FP  -  0.10,  RTn  —  40:  then  from  Table  2,  CEP  =  23.717  <  36; 
can  we  conclude  the  contractor  should  be  rewarded?  The  answer  is  “yes”  because,  from 
Table4.  36  is  to  the  right  of  the  interval  (18.3685.  33.4977),  equivalently  greater  than  CHR. 
which  in  90%  of  similar  test  cases  will  contain  the  CEPT  and  will  be  less  than  36.  Instead, 
suppose  the  RTn  -  16.  Since  the  CEP  =  23.717  is  greater,  we  try  to  judge  if  the  contractor 
should  be  penalized.  Namely,  is  16*  1.1  =  17.6  <  CLR(=  18.3685)?  The  answer  is  again 
“yes”  in  the  numerical  example,  since  in  90%  of  repeated  like  tests  the  CEPT  would  exceed 
17.6.  In  either  case,  if  the  answer  was  negative  due  to  a  different  choice  of  RTn,  then  more 
testing  would  be  indicated  (MOT). 

Before  considering  the  procedure  of  more  testing,  it  is  important  to  arrive  at  a  conclusion 
with  a  minimum  of  testing.  Hence  we  ask  in  the  COR  case  if  we  can  lower  CHR;  i.e.,  ease 
the  requirement  for  a  decision.  Therefore,  using  DGINV,  we  ask  for  the  value  of  \2  such 


C(Xa,  oo,  f)  =  1  -  a  =  0.9486833,  C(0,  f)  =  a  =  0.0513167.  (34) 

so  that. 

f  s l/al  >  xl  =>  the  Cl  for  <rx  :  0  <  <rx  <  b  sx,  <ry  >0,  b  =  y/f Txl  >  1.  (35) 

f  Sy/o'y  >  xi  =>  the  Cl  for  ay  :  0  <  oy  <  b  sy,  ax  >  0,  b  =  y/f/Xa  >  1-  (36) 

The  intersection  of  the  two  infinite  CIs  is  shown  as  the  hatched  region  in  Figure  5. 

Note  then,  as  in  Section  IV,  the  (1  -  a)2  Cl  for  the  CEPT  is  given  by  solving  for  CHR 

K  =  CEP/sv  =  CHR/fbsy)  =►  CHR  =  b  CEP.  (37) 

where,  in  general, 

CHR  <  CHR  (see  Section  IV  for  CHR).  (38) 

In  the  case  where  (2)  holds,  INVELP  with  standard  deviations  bsx  and  bsy  is  used  to  find 


CHR.  For  the  same  reasons  given  in  the  previous  section  in  the  discussion  of  the  nonzero 
means  case,  (37)  cannot  be  used  (See  the  paragraph  before  Table 4.). 

We  ask  in  the  COP  case  if  we  can  raise  CLR;  i.e.,  ease  the  requirement  for  a  decision. 
Therefore,  using  DGINV,  we  ask  for  the  value  of  xfi~a)  su°h  that 


that 

C(0,  Xq-o) 

,  f)  =  1  -  Q  =  0.9486833, 

(39) 

f  <  xl-a) 

=>  the  Cl  for  ax  : 

0  <  a  sx  <  ax,  ay  >  0,  a  =  Vf/xfi-a)  <  1- 

(40) 

f  4/av  ^  X(i-Q) 

=>  the  Cl  for  ay  : 

0  <  a  sy  <  cry,  (tx  >  0,  a=\/f/x?i _Q)  <  L 

(41) 
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The  intersection  of  the  two  infinite  CIs  is  shown  as  the  hatched  region  in  Figure  6. 

Note  then,  as  in  Section  IV,  the  (1  -  a)2  Cl  for  the  CEPT  is  given  by  solving  for  CLR 

K  =  CEP/sx  =  CLR/(asx)  =►  CLR  =  a  CEP,  (42) 

whore,  in  general,  CLR  <  ^  (gee  Section  ,v  for  CLR).  (43) 

In  the  case  where  (2)  is  used  instead  of  (1),  INVELP  with  standard  deviations  asx  and 
asy  is  used  to  find  CLR.  Equation  (42)  cannot  be  used  for  the  same  reasons  given  in  the 
previous  section;  (see  the  paragraph  before  Table  4.) 


Figure 5.  With  Confidence  (1  —  a)2,  0  <  RT  <  bR 


Figure6.  With  Confidence  (1  —  a)2,  aR  <  RT  <  oo 
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Table  5  ContainsCLR  And  CHR;  It  Supplements  Table 4. 


Table  5.  Listing  of  CLR  >  CLR.  CHR  <  CHR 
Using  Data  From  Tables  1  and  4 


CI|  Means  -► 

x  =  0.  y  =  0 

x  ^  0,  y^O 

CLR 

CHR 

19.1306 

31.5975 

19.5547 

32.5010 

Perhaps  the  most  efficient  strategy  in  what  has  been  discussed  so  far,  to  evaluate  whether 
COR,  COP.  or  MOT  holds,  is  to  determine  CLR  or  CHR  at  each  new  set  of  sequential 
bombing  tests.  Given  RTn,  FP.  n().  An.  and  Jmax,  where  nO  denotes  the  initial  number  of 
bomb  drops.  An  is  the  increment  number  of  bomb  drops  added  at  each  test  after  the  initial 
nO  drops,  and  Jmax  denotes  the  assigned  number  of  Ans,  incremental  bomb  drops  to  be 
made,  after  which  all  testing  stops. 

Using  the  data  from  Tables  1  and  2,  with  (1).  CEP  =  23.717  and  with  (2),  CEP  =  24.209, 
assign  the  following: 

nO  =  10,  An  =  4,  Jmax  =  3.  (44) 

Tests,  as  described  in  Section  4  and  this  section,  are  conducted  for  each  n,  where 


n  =  nO  +  j An,  j  =  0,  1,  • •  •  ,  Jmax. 


(45) 


The  results  are  given  in  Table  6  of  page  15.  CIs  can  be  identified  from  the  table  accordingly: 


1.  Conventional  (1  —  a)2  (=  .90)  Cl  for  the  CEPT 

2.  Minimum  (1  —  a)2  (=  .90)  Cl  for  the  CEPT 

3.  Symmetric  (1  —  a)2  (=  .90)  Cl  for  the  CEPT 

4.  Maximum  CLR  (1  -  a)2  (=  .90)  Cl  for  the  CEPT 

5.  Minimum  CHR  (1  -  a)2  (=  .90)  Cl  for  the  CEPT 


(L(l),  H(l))  =  (CLR.  CHR) 
(L(2),  H(2))  =  (MLR.  MHR) 
(L(3),  H(3))  =  (SLR,  SHR) 
(L(4),  H(4))  =  (CLR,  oo ) 
(L(5),  H(5))  =  (0,  CHR) 


For  example,  for  the  zero  means  case,  for  the  symmetric  Cl  at  n  =  18,  one  obtains  from 
the  table  (13.7266.  30.6586)  with  a  CEP  =  22.1926;  the  minimum  CHR  Cl  is  (0,  30.6462). 
For  the  nonzero  means  case,  one  obtains  for  the  symmetric  Cl  at  n  =  18,  (14.0863.  31.6683) 
with  a  CEP  =  22.8205;  the  minimum  CHR  =  Cl  is  (0,  31.6598). 

If  RTn  =  16.5  and  FP  =  0.10,  then  COP  holds  for  the  zero  means  case,  since 
18.15  <  CLR  =  19.1306  <  CEP  at  n  =  22.  In  the  nonzero  means  case,  COP  also  holds 
since  18.15  <  CLR  =  18.8265  <  CEP  for  n  =  10.  On  the  other  hand,  let  RTn  =  34.2  and 
FP  =  0.1,  then  COR  holds  in  the  zero  means  case,  since  CHR  =  30.6462  <  30.78  for  n  =  18, 
but  in  the  nonzero  means  cases  MOT  holds,  since  for  n  =  10,  14,  18,  22,  CHR  >  30.78. 

In  the  unlikely  event  that  RP  «  CEPT  or  RN  «  CEPT,  it  may  require  a  very  large  value 
of  n  to  conclude  whether  COP  or  COR  holds. 
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Table  6.  Listing  of  CIs  for  CEPT  Using  Data  From  Previous  Tables 


L(1)=CLR, 

L(2)=MLR,  L(3)=SLR, 

L(4)=CLR,  L(5)= 

4) 

H(1)=CHR, 

H(2)=MHR,  H(3)=SHR, 

H(4)=oo, 

H(5)=CHR 

j,  n,  m 

L(m) 

H(m) 

CEP 

j,  n,  m 

L(m) 

H(m) 

CEP 

X 

y 

Sx 

sy 

X 

y 

Sx 

sy 

0.0 

0.0 

29.4891 

12.4880 

3.8299 

-2.1856 

30.8210 

12.9603 

0,  10,  1 

16.9769 

42.4081 

24.2501 

0,  10,  1 

17.7882 

46.1287 

25.5357 

0,  10,  2 

15.6089 

39.3400 

24.2501 

0,  10,  2 

16.2745 

42.4788 

25.5357 

0,  10,  3 

10.0136 

38.4866 

24.2501 

0,  10,  3 

9.8215 

41.5872 

25.5357 

0.  10,  4 

17.9640 

00 

24.2501 

0,  10,  4 

18.8265 

00 

25.5357 

0,  10,  5 

0.0 

38.4866 

24.2501 

0,  10,  5 

0.0 

41.5872 

25.5357 

X 

y 

sx 

sy 

X 

y 

Sx 

sy 

0.0 

0.0 

26.9855 

13.4646 

1.5335 

-2.3071 

27.9590 

13.7662 

1,  14,  1 

17.2129 

36.9126 

23.4708 

1,  14,  1 

17.7751 

38.9652 

24.3463 

1,  14,  2 

16.2033 

34.9524 

23.4708 

1,  14,  2 

16.6815 

36.7593 

24.3463 

1,  14,3 

12.7839 

34.1578 

23.4708 

1,14,3 

12.8455 

35.9331 

24.3463 

1,14,4 

18.0811 

OC 

23.4708 

1,14,4 

18.6803 

DO 

24.3463 

1,14,5 

0.0 

34.1564 

23.4708 

1,  14,  5 

0.0 

35.9326 

24.3463 

X 

y 

Sx 

sy 

X 

y 

Sx 

Sy 

0.0 

0.0 

26.2413 

12.0934 

4.0512 

-2.2699 

26.6783 

12.2229 

2,  18,  1 

16.7947 

32.7412 

22.1926 

2,  18,  1 

17.3269 

33.8762 

22.8205 

2,  18,  2 

16.0203 

31.3648 

22.1926 

2,  18,  2 

16.5234 

32.3927 

22.8205 

2,  18,  3 

13.7266 

30.6586 

22.1926 

2,  18,  3 

14.0863 

31.6683 

22.8205 

2,  18,  4 

17.5552 

OC 

22.1926 

2,  18,  4 

18.0945 

OC 

22.8205 

2,  18,  5 

0.0 

30.6462 

22.1926 

2,  18,  5 

0.0 

31.6598 

22.8205 

X 

y 

sx 

sy 

X 

y 

Sx 

Sy 

0.0 

0.0 

27.5922 

13.3219 

0.2657 

-2.4986 

28.2403 

13.3934 

3,  22,  1 

18.3685 

33.4977 

23.7170 

3,  22,  1 

18.7754 

34.5007 

24.2689 

3,  22,  2 

17.6720 

32.3322 

23.7170 

3,  22,  2 

18.0454 

33.2535 

24.2689 

3,  22,  3 

15.7998 

31.6342 

23.7170 

3,  22,  3 

16.0451 

32.5316 

24.2689 

3,  22,  4 

19.1306 

OC 

23.7170 

3,  22,  4 

19.5547 

OO 

24.2689 

3,  22,  5 

0.0 

31.5975 

23.7170 

3,  22,  5 

0.0 

32.5010 

24.2689 

The  next  section  describes  the  Fortran  95  computer  program  that  outputs  the  results  of 
Table  6  as  well  as  whether  MOT  or  COP  or  COR  holds  given  RN,  RP. 
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VI.  FORTRAN  PROGRAMS  FOR  CEP  AND  CONFIDENCE  INTERVALS 

The  total  software  package  is  made  up  of  69  double-precision  Fortran  95  subprograms. 
Each  of  their  call  lines  is  listed  in  Appendix  B.  The  first  14  were  developed  by  the  author;  the 
remaining  55  are  supporting  routines  from  the  NSWC  Library  of  Mathematics  Subroutines. 

[6].  Of  the  first  14.  we  shall  discuss  the  first  5;  the  remaining  9  are  supporting  routines.  The 
software  is  stored  in  a  file  named  WEEV4.FOR. 

The  Program  routine,  MEANSTDV,  is  actually  not  necessary  and  can  be  deleted.  It  is  a 
test  routine  for  the  software  that  uses,  as  input,  a  file  JDATA.  If  the  user  wishes,  the  software 
could  be  run  with  a  read-in  input  of  0.5  or  0.9.  and  2  when  requested,  and  directing  the 
results  by  WEEV4  >  f8;  then  making  a  DOS  file  “compare,”  FC  f8  JDATATEST.  The  files 
should  be  duplicates  except  for  the  last  line,  which  calls  on  the  Lahey  subroutines  DATE 
and  TIME.  The  required  JDATA  on  the  first  line  would  list  the  two  integers  1  and  22,  on  the 
2nd  line  store  xi,  yj,  and  the  following  21  lines  would  contain  the  remaining  miss  distances 
of  Table  1,  x,,  yj,  (i  =  2,  22). 

The  2nd  listed  routine, 

CONFREG(N,X.Y.PCEP,PC,ISV,RN,RP,XYM,SXY,RCEP.RLH.DEC), 
is  the  main  routine  for  WEEV4,  and  it  is  the  only  one  a  user  needs  to  call.  Its  call  line 

specifies  the  necessary  input  and  output.  The  first  ten  variables  are  required  input;  the  last 

three  are  output. 

INPUT 

N  Number  of  bomb  drops. 

X  N-element  array  containing  the  x-coordinate  miss  distances. 

Y  N-element  array  containing  the  y-eoordinate  miss  distances. 

PCEP  Probability  measure  for  the  CEPT.  Usually  taken  as  0.50. 

PC  Assigned  confidence  level  for  CIs.  Taken  as  0.9  in  the  examples. 

ISV  If  0(>  0),  then  do  zero(nonzero)  means  case. 

RN  Value  (1  —  FP)  RTn.  In  computer  numerical  example  FP  =  0.125.  RTn  =  162/3, 
RN  —  14.583333.  If  RN  <  0,  the  contractor  evaluation  is  skipped. 

RP  Value  (1  +  FP)  RTn.  In  computer  example  RP  =  18.75. 

XYM  Two-element  array  containing  means  x  and  y. 

SXY  Two-element  array  containing  standard  deviations  sx  and  sy. 
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OUTPUT 

RCEP  Estimated  CEP.  If  ISV  =  0  uses  INVGCE,  else  uses  INVELP. 

RLH  Eight-element  array  with  the  following  elements  (See  Table 6): 

RLH(l)  —  CLR,  Left  endpoint  of  the  conventional  Cl  for  CEPT. 

RLH (2)  -  CHR,  Right  endpoint  of  the  conventional  Cl  for  CEPT. 

RLH (3)  —  MLR.  Left  endpoint  of  the  minimum  Cl  for  CEPT. 

RLH (4)  —  MHR.  Right  endpoint  of  the  minimum  Cl  for  CEPT. 

RLH (5)  —  SLR,  Left  endpoint  of  the  symmetric  Cl  for  CEPT. 

RLH(6)  —  SHR.  Right  endpoint  of  the  symmetric  Cl  for  CEPT. 

If  CHR  >  2,  then  SLR  and  SHR  set  to  -1. 

RLH(7)  —  CLR,  Maximum  CLR,  left  endpoint  of  a  Cl  for  CEPT. 

RLH(8)  —  CHR,  Minimum  CHR,  right  endpoint  of  a  Cl  for  CEPT. 

DEC  A  3-character  variable  specifying  contractor  performance: 

COR.  COP,  MOT;  if  input  RN  <  0,  ignore  DEC. 

The  next  three  routines  are  called  by  CONFREG.  MNSTD  is  used  to  compute  x,  y,  and 
sx,  sy.  CONFMS  is  used  to  find  the  three  CIs:  namely,  CCR,  MCR,  and  SCR.  EVAL  is  used 
to  determine  if  the  contractor  is  to  be  rewarded,  COR;  if  the  contractor  is  to  be  penalized, 
COP;  or  if  more  testing  is  needed,  MOT. 
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APPENDIX  A 


GLOSSARY 

The  first  number  listed  with  each  item,  refers  to  the  page  where  it  is  first  introduced. 


C(7,  <5,  f) 

CHR 

CLR 

CEPT 

a 

An 

(1-q) 

CCI 

CCR 

CEP 

CHR 

CHx 

CHy 

Cl 

CLR 

CLx 

CLy 

COP 

COR 

ELP 

ELP 

f 

FP 

GCE 

GCE 

INVELP 

INVGCE 

Jrriax 

K 

K 


Integral  of  the  Chi-square  distribution  density  from  7  to  S,  7 
The  minimum  value  CHR  can  take  for  a  given  a,  12 
The  maximum  value  CLR  can  take  for  a  given  a,  13 
True  CEP,  1,  6 
Confidence  level,  7,  9 

Increment  number  of  bomb  drops  added  at  each  test  after  the  first,  14 

(1  -  a)  =  VPC,  7 

Conventional  confidence  interval,  6 

Conventional  confidence  interval  for  CEPT,  11 

Estimated  circular  error  probable,  1 

R,  Right  endpoint  of  the  CCI,  CCR,  for  CEPT,  11 

Right  endpoint  of  the  CCI  for  <rx.  See  (24),  7 

Right  endpoint  of  the  CCI  for  ay.  See  (24),  7 

Confidence  interval (s),  1,  8 

R,  Left  endpoint  of  the  CCI,  CCR,  for  CEPT,  11 

Left  endpoint  of  the  CCI  for  crx.  See  (24),  7 

Left  endpoint  of  the  CCI  for  ay.  See  (24),  7 

Contractor  is  penalized,  12 

Contractor  is  rewarded,  12 

Elliptic  Coverage  Function,  2 

ELP  subroutine  to  compute  ELP,  5 

number  of  degrees  of  freedom,  2 

A  fixed  decimal  associated  with  RTn,  11 

GCE  subroutine  to  evaluate  GCE,  3 

Generalized  Circular  Error  Function,  2 

Subroutine  to  compute  R  of  P(R,  H,  K,  sx,  sy),  5 

Subroutine  to  compute  K  of  V(K,c),  4 

Maximum  number  of  Ans,  14 

Dimensionless  radius  of  GCE.  K  =  R/u,  3 

y-coordinate  of  the  mean  of  the  Elliptic  Coverage  Function,  5 
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MCI 

- 

Minimized  confidence  interval,  G.9 

MCR 

- 

Minimum  confidence  interval  for  CEPT,  1 1 

MHR 

- 

R,  Right  endpoint  of  the  MCI.  MCR,  for  CEPT,  11 

MHx 

- 

Right  endpoint  of  the  MCI  for  ax.  See  (33),  10 

MHy 

- 

Right  endpoint  of  the  MCI  for  ay.  See  (33),  10 

MLR 

- 

R,  Left  endpoint  of  the  MCI,  MCR  .  for  CEPT,  1 1 

MLx 

- 

Left  endpoint  of  the  MCI  for  ax.  See  (33),  10 

MLy 

- 

Left  endpoint  of  the  MCI  for  ay.  See  (33),  10 

MOT 

- 

More  testing  required,  12 

nO 

- 

Initial  number  of  bomb  drops,  14 

PC 

- 

Assigned  probability  for  the  confidence  interval  of  the  CEPT.  7 

PCEP 

- 

Percentage/ 100  of  the  UBND  contained  in  the  target  circle,  1 

P(R.  H.  K,  sx,  sy) 

- 

Elliptic  Coverage  Function,  2 

R 

- 

Estimated  CEP  (computed),  1 

R 

- 

Used  interchangeably  with  the  estimated  CEP,  2 

RED 

- 

Radar  bomb  director,  1 

RN 

- 

=(l-FP)RTn.  12 

RP 

- 

=  (l+FP)RTn,  12 

RT 

- 

True  CEP=CEPT,  1 

RTn 

- 

A  pre  specified  nominal  value  for  the  CEP,  1 

SCI 

- 

Symmetrical  confidence  interval,  6, 9 

SCR 

- 

Symmetrical  confidence  interval  for  CEPT,  11 

SHR 

- 

R,  Right  endpoint  of  the  SCI.  SCR.  for  CEPT,  11 

SHx 

- 

Right  endpoint  of  the  SCI  for  <rx.  See  (33),  10 

SHy 

- 

Right  endpoint  of  the  SCI  for  rrv.  See  (33),  10 

SLR 

- 

R,  Left  endpoint  of  the  SCI,  SCR.  for  CEPT,  1 1 

SLx 

- 

Left  endpoint  of  the  SCI  for  (tx.  See  (33),  10 

SLy 

- 

Left  endpoint  of  the  SCI  for  ay.  See  (33),  10 

UBND 

- 

Uncorrelated  Bivariate  Normal  Distribution,  1 

V(K.c) 

- 

Generalized  Circular  Error  Function,  2 
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APPENDIX  B 

LIST  OF  FORTRAN  95  SUBPROGRAMS 

1  PROGRAM  MEANSTDV 

2  SUBROUTINE  CONFREG(N,X,Y,PCEP,PC,ISV,RN,RP,XYM,SXY,RCEP,RLH,DEC) 

3  SUBROUTINE  MNSTD(N,X,XM,VX2,SX2) 

4  SUBROUTINE  CONFMS(N,P.RLO.RHO,RL1,RH1,RL2.RH2) 

5  SUBROUTINE  EVAL(N,RN,RP,CHR) 

6  DOUBLE  PRECISION  FUNCTION  FD(X)  IFOR  MINIMUM  CONF  INTERVALS, 

USED  IN  FMIN  IN  CONFMS 

7  DOUBLE  PRECISION  FUNCTION  FA(X)  IFOR  SYMMETRIC  CONF  INTERVALS, 

USED  IN  FMIN  IN  CONFMS 

8  SUBROUTINE  INVGCE(P,C,R,IJ,IERR) 

9  SUBROUTINE  INVELP(PP,H1,K1,S1,S2,R,P,I) 

10  DOUBLE  PRECISION  FUNCTION  F2(X)  INEEDED  FOR  DZERO  OF  INVELP 

11  SUBROUTINE  DQ(R,PX,II) 

12  DOUBLE  PRECISION  FUNCTION  F1(X)  IFUNCTION  FOR  DQXGS  ROUTINE 

13  SUBROUTINE  GRUB(R)  IINITIAL  ESTIMATE  FOR  R  FROM  GRUBBS’  APPROX 

14  SUBROUTINE  PSQR(A.P,X,I)  IlST  APPROX,  X,  USES  GRUBBS’  ESTIMATE 

15  INTEGER  FUNCTION  IPMPAR  (I) 

16  DOUBLE  PRECISION  FUNCTION  DPMPAR  (I) 

17  DOUBLE  PRECISION  FUNCTION  DEPSLN  (L) 

18  DOUBLE  PRECISION  FUNCTION  DXPARG  (L) 

19  DOUBLE  PRECISION  FUNCTION  DSIN1  (X) 

20  DOUBLE  PRECISION  FUNCTION  REXP  (X) 

21  DOUBLE  PRECISION  FUNCTION  DREXP  (X) 

22  DOUBLE  PRECISION  FUNCTION  ALNREL(A) 

23  DOUBLE  PRECISION  FUNCTION  DLNREL  (A) 

24  DOUBLE  PRECISION  FUNCTION  RLOG(X) 

25  DOUBLE  PRECISION  FUNCTION  DRLOG  (X) 

26  DOUBLE  PRECISION  FUNCTION  ERF  (X) 

27  DOUBLE  PRECISION  FUNCTION  ERFCl  (IND,  X) 

28  DOUBLE  PRECISION  FUNCTION  DERF  (X) 

29  DOUBLE  PRECISION  FUNCTION  DERFC  (X) 

30  DOUBLE  PRECISION  FUNCTION  DERFC  1  (IND,  X) 

31  DOUBLE  PRECISION  FUNCTION  DERFCO  (X) 

32  DOUBLE  PRECISION  FUNCTION  ERFI  (P,  Q) 
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33  DOUBLE  PRECISION  FUNCTION  DERFI  (P.  Q) 

34  DOUBLE  PRECISION  FUNCTION  DAERF  (X,  H) 

35  SUBROUTINE  PNI  (P.  Q.  D.  W,  IERR) 

36  SUBROUTINE  DPNI  (P.  Q.  D.  W.  IERR) 

37  DOUBLE  PRECISION  FUNCTION  GAMMA(A) 

38  DOUBLE  PRECISION  FUNCTION  GLOG(X) 

39  DOUBLE  PRECISION  FUNCTION  GAM1(A) 

40  DOUBLE  PRECISION  FUNCTION  GAMLN  (A) 

41  DOUBLE  PRECISION  FUNCTION  GAMLN  1  (A) 

42  DOUBLE  PRECISION  FUNCTION  DGAMMA(A) 

43  DOUBLE  PRECISION  FUNCTION  DPDEL(X) 

44  DOUBLE  PRECISION  FUNCTION  DGAM1  (X) 

45  DOUBLE  PRECISION  FUNCTION  DGAMLN  (A) 

46  DOUBLE  PRECISION  FUNCTION  DGMLN1  (X) 

47  SUBROUTINE  G RAT  10  (A,  X,  ANS,  QANS,  IND) 

48  DOUBLE  PRECISION  FUNCTION  RCOMP  (A,  X) 

49  DOUBLE  PRECISION  FUNCTION  DRCOMP  (A.  X) 

50  SUBROUTINE  GAMINV  (A,  X,  XO,  P,  Q,  IERR) 

51  SUBROUTINE  DGRAT  (A.  X,  ANS,  QANS,  IERR) 

52  SUBROUTINE  DGR29  (A,  Y,  L,  Z,  RTA,  ANS,  QANS) 

53  SUBROUTINE  DGR17  (A,  Y,  L,  Z,  RTA,  ANS,  QANS) 

54  SUBROUTINE  DGINV  (A,  X,  P,  Q,  IERR) 

55  SUBROUTINE  BESI(X.  ALPHA,  KODE,  N,  Y,  NZ) 

56  SUBROUTINE  ASIK(X,FNU,KODE,FLGIK,RA.ARG,IN,TOL,Y) 

57  SUBROUTINE  CIRCV  (R,  D.  J,  P,  IERR) 

58  SUBROUTINE  ERFCO  (IND.  X,  E,  Y) 

59  DOUBLE  PRECISION  FUNCTION  ERFCR  (X) 

60  DOUBLE  PRECISION  FUNCTION  DZERO  (F,  AX,  BX,  AERR,  RERR) 

61  SUBROUTINE  FMIN(F,  AO,  BO,  X,  W,  AERR,  RERR,  ERROR,  IND) 

62  SUBROUTINE  DQPSRT(LIMIT,LAST,MAXERR.ERMAX,ELIST.IORD,NRMAX) 

63  SUBROUTINE  DQELG  (N,  EPSTAB,  RESULT,  ABSERR,  RES3LA,  NRES, 

EPMACH.  OFLOW) 

64  SUBROUTINE  DQXGS  (F,A.B,EPSABS.EPSREL, RESULT, ABSERR.IER. 

LIMIT.LENIW,LENW, LAST, IWORK, WORK) 

65  SUBROUTINE  DQXGSE(F, A, B,EPSABS,EPSREL,LIMIT, RESULT, ABSERR, 

IER,ALIST,BLIST,RLIST,ELIST,IORD,LAST,VALP,VALN,LP,LN) 


B-4 


NSWCDD/TR— 07/13 


GO  SUBROUTINE  DQXCPY  (A,  B,  L) 

67  SUBROUTINE  DQXLQM  (F,  A,B, RESULT,  ABSERR,RESABS.RESASC,VR,  VS, LR  ,LS, 

KEY,EPMACH,UFLOW,OFLOW) 

68  SUBROUTINE  DQXRUL  (F,XL.XU,Y.YA,YM,KE,K1.FV1.FV2,L1,L2) 

69  SUBROUTINE  DQXRRD  (F,Z,LZ,XL,XU,R,S,LR,LS) 
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